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Abstract -We present a continuum model of the coupling between cells and substrate that 
accounts for some of the observed substrate-stiffness dependence of cell properties. The cell 
is modeled as an elastic active gel, adapting recently developed continuum theories of active 
viscoelastic fluids. The coupling to the substrate enters as a boundary condition that relates the 
cell's deformation field to local stress gradients. In the presence of activity, the coupling to the 
substrate yields spatially inhomogeneous contractile stresses and deformations in the cell and can 
enhance polarization, breaking the cell's front-rear symmetry. 
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• Introduction. — Many cell properties, including cell 
. shape, migration and differentiation, are critically con- 
trolled by the strength and nature of the cell's adhesion to 
Qpa solid substrate and by the substrate's mechanical prop- 
'— erties [lj. For instance, it has been demonstrated that 
cell differentiation is optimized in a narrow range of ma- 
trix rigidity [2] and that the stiffness of the substrate can 
direct lineage specification of human mesenchymal stem 
O I cells [3j. In endothelial cells, adhesion to a substrate plays 
a crucial role in guiding cell migration and controlling a 
^-".number of physiological processes, including vascular de- 
velopment, wound healing, and tumor spreading [4j. Fi- 
broblasts and endothelial cells seem to generate more trac- 
l tion force and develop a broader and flatter morphology 
■ • on stiff substrates than they do on soft but equally adhe- 
^ sive surfaces [5j[6] . They show an abrupt change in their 
spread area within a narrow range of substrate stiffnesses. 
1— j This spreading also coincides with the appearance of stress 
fibers in the cytoskeleton, corresponding to the onset of a 
substantial amount of polarization within the cell [6] . Fi- 
nally, such cells preferentially move from a soft to a hard 
surface and migrate faster on stiffcr substrates [7j. The 
mechanical interaction of cells with a surrounding matrix 
is to a great extent controlled by contractile forces gener- 
ated by interactions between filamentary actin and myosin 
proteins in the cytoskeleton. Such forces are then trans- 
mitted by cells to their surroundings through the action 
of focal adhesions that produce elastic stresses both in the 
cell and in the surrounding matrix. Cells in turn are ca- 
pable of responding to the substrate stiffness by adjusting 



their own adhesion and elastic properties, with important 
implications for cell motility and shape 1 1 , 



In this letter we present a simple model of the coupling 
between cells and substrate that accounts for some of the 
observed substrate-stiffness dependence of cell properties. 
The cell itself is modeled as an elastic active gel, adapting 
recently developed continuum theories of active viscoelas- 
tic fluids [9 11 . In these models the transduction of chem- 
ical energy from ATP hydrolysis into mechanical work by 
myosin motor proteins pulling on actin filaments yields 
active contractile contributions to the local stresses. The 
continuum theory of such active liquids has led to several 
predictions, including the onset of spontaneous deforma- 
tion and flow in active films 12p3 and the retrograde flow 
of actin in the lamellipodium of crawling cells [TT] . Active 
liquids cannot, however, support elastic stresses at long 
times, as required for the understanding of the crawling 
dynamics of the lamellipodium and of active contractions 
in living cells. Models of active elastic solids on the other 
hand have been shown to account for the contractility 
and stiffening of in- vitro actomyosin networks 14f(16| and 
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the spontaneous oscillations of muscle sarcomeres 
Very recently a continuum model of a one-dimensional po- 
lar, active elastic solid has also been used to describe the 



alternating polarity patterns observed in stress fibers 19 



In all these cases the elastic nature of the network at low 
frequency is crucial to provide the restoring forces needed 
to support deformations and oscillatory behavior. 

We model a cell as an elastic active film anchored to a 
solid substrate and study the static response of the film 
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to variations in the strength of the anchoring. Although 
in the following we refer to our system as a cell, we stress 
that, on different length scales, the active elastic gel could 
also serve as a model for a confluent cell monolayer on a 
substrate. The coupling of the cell to the substrate enters 
via a boundary condition controlled by a "stiffness" pa- 
rameter that depends on both the cell/substrate adhesion 
as well as the substrate rigidity. The description is macro- 
scopic and applies on length scales large compared to the 
typical mesh size of the actin network in the cell lamel- 
lipodium (or large compared to the typical cell size in the 
case of a cell monolayer). By solving the elasticity and 
force balance equations in a simple one-dimensional ge- 
ometry we obtain several experimentally relevant results. 
First, in an isotropic active gel substrate anchoring yields 
stresses and contractile deformations. The stress and de- 
formation profiles for an isotropic active elastic gel are 
shown in the top frame of Fig. [I] The stress is largest at 
the center of the cell. Interestingly, a very similar profile 
of tensile stresses has been observed in confluent mono- 
layers of migrating epithelial cells 20 , where the stress 



increases as a function of the distance from the leading 
edge of the migrating layer and reaches its maximum at 
the center of the cell colony. Although our model considers 
stationary active elastic layers (and the resulting stresses 
are contractile as opposed to tensile), in both cases these 
stresses originate from active processes in the cell, driven 
by ATP consumption. The deformation of the active layer 
is largest at the cell boundaries (see Fig [I] top frame), 
as seen in experiments imaging traction forces exerted by 
cells on substrates [2l] and its overall magnitude increases 
with cell activity. The density of the active gel layer is 
concentrated at the boundary, where the local contractile 
deformations are largest. The net deformation of the cell 
over its length is shown in the bottom frame of Fig. [I] and 
it increases monotonically with decreasing substrate stiff- 
ness, in qualitative agreement with experiments on fibrob- 
lasts showing that these cells are more extended on stiff 
substrates |6]. Finally, if the cell is polarized on average, 
the coupling to the substrate generates a spatially inho- 
mogeneous polarization profile inside the cell. The mean 
polarization is enhanced over its value in the absence of 
substrate anchoring and it is a non-monotonic function of 
substrate stiffness (see Fig. |4| . This result is in qualita- 
tive agreement with recent experiments that have demon- 
strated an intimate relation between the matrix rigidity 
and the alignment of cell fibers within the cell, suggesting 
that maximum alignment may be obtained for an optimal 



value of the substrate rigidity 22 



The active gel model. — The cell is modeled as an 
active gel described in terms of a density, p(r,t), and a 
displacement field, u(r, t), characterizing local deforma- 
tions. In addition, to account for the possibility of cell 
polarization as may be induced by directed myosin motion 
and/or filament treadmilling, we introduce a polar orien- 
tational order parameter field, P(r,i). Although we are 
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Fig. 1: Top: stress a(x)/(A/j, (dashed line) and deformation 
u(x)B /(,A[i (solid line) profiles as functions of the position x 
inside a cell of length L for X/L — 0.25. Bottom: the cell's 
total deformation Al = u(0) — u(L) as a function of X/L. In 
the plot the deformation Al is normalized to its maximum 
value (A/i/B. 



describing a system out of equilibrium, it is convenient to 
formulate the model in terms of a local free energy density 
/ = fel + fp + fw, with 
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Here f e i is the energy of elastic deformations, with B and 
G the compressional and shear elastic moduli of the gel, 
respectively, u 
tensor, with Un — 



= ^{diUj + djUi) the symmetrized strain 
1 6aUkk and d the dimensionality. 



The first two terms in Eq. (Ib), with b > 0, allow the on- 



set of a homogeneous polarized state when a < 0; the last 
term is the energy cost for spatially inhomogeneous de- 
formations of the polarization. We have used an isotropic 
elastic constant approximation, with K a stiffness param- 
eter characterizing the cost of both splay and bend defor- 
mations. Finally, the contribution f w couples strain and 
polarization and is unique to polar systems 13 19. It 



describes the fact that in the active polar system consid- 
ered here, like in liquid crystal elastomers, a local strain 
is always associated with a local gradient in polarization. 
Such gradients will align or oppose each other depending 
on the sign of the phenomenological parameters w and w' , 
which are controlled by microscopic physics. A positive 
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sign indicates that an increase of density is accompanied 
by positive splay (or enhanced polarization in one dimen- 
sion). In active actomyosin systems filament polarity can 
be induced by both myosin motion and by treadmilling. 
If the polarization is defined as positive when pointing to- 
wards the plus (barbed) end of the filament, i.e., the direc- 
tion towards which myosin proteins walk, the forces trans- 
mitted by myosin procession will yield filament motion in 
the direction of negative polarization, corresponding to 
w < 23 . In contrast, treadmilling, where polarization 
occurs at the barbed end, corresponds to w > 0. Density 
variations Sp = p — pq from the equilibrium value, po, are 
slaved to the local strain according to Sp/po = — V • u. 
The stress tensor is written as the sum of reversible and 



active contributions as cry = cr|"- + 
The two contributions are given by 
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of, = C(p)AM (J + C a Ap PiPj 



Active stresses arise because the gel is driven out of equi- 
librium by continuous input of energy from the hydrolysis 
of ATP, characterized by the chemical potential difference 
Ap between ATP and its products. For simplicity, we as- 
sume here A/z to be constant, although situations where 
inhomogeneities in Ap may arise, for instance, from in- 
homogeneous myosin distribution within the actin lamel- 
lipodium are also of interest. The experimentally observed 
contractile effect of myosin corresponds to positive values 
of the coefficients £ and ( a , that characterize the isotropic 
and anisotropic stress per unit A/i, respectively, due to 
the action of active myosin crosslinkers 10 24 25 . In 



polar gels there are also active stresses proportional to 
Api^P, 



djPi) 
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We neglect these terms here as 
terms of similar structure already arise from the coupling 
terms in f w . By letting p = p — p "V ■ u, we can write 
((p)Ap ~ CoAp — QiApukk- The second term describes 
active renormalization of the compressional modulus B 
of the gel and can yield a contractile instability 17 18]. 



These effects have been described elsewhere 18 and will 
not be discussed here, where we will assume we are in 
a regime where the gel is elastically stable. Finally, we 
note that the parameters a, w and w' may also in general 
depend on Ap as cell polarity is induced by ATP-driven 
processes. For simplicity we keep these parameters fixed 
below. 

Force balance requires 



djO-ij 







(3) 



The coupling to the substrate (assumed for simplicity 
isotropic) is introduced as a boundary condition 27 by 



requiring [uy^L = Eui(r = r s ), where n is a unit nor- 
mal to the substrate and both sides of the equation are 
evaluated at points r s on the substrate. Although in the 
following we will often refer to the parameter E as the 



substrate stiffness, it is important to keep in mind that 
E is controlled not only by the substrate rigidity, but 



also by the properties of the cell/substrate adhesions 28 



Anisotropic substrates are not considered here, but can be 
described by a generalized boundary condition where E is 
a tensor quantity and will be discussed in a later publica- 
tion. Finally, variations in the polarization are described 
by the equation 



d t P + /3 (P • V) P = Th 



(4) 



with j3 an advective coupling arising from ATP driven pro- 
cesses, such as treadmilling [13||26| , T an inverse friction, 
and h = — |^ the molecular field, given by 



■b\P\ 2 )P l + KV 2 P l 



wdjiitj + w'diu kk . (5) 



Here f3 is an active velocity and is controlled by the activity 
Ap. In the following we write j3/(LT) = (pAp, with L the 
typical size of the active gel. 

Isotropic cell. — We begin by considering the case of 
an isotropic cell and neglect the coupling to polarization. 
For simplicity, we consider a quasi-one-dimensional model 
where the cell is a thin sheet of active gel of thickness 
z extending from x — to x — L, with L » h. The 
substrate is flat and located at z = 0. Although this is 
of course a gross simplification, we will see below that it 
captures the substrate-induced stresses and deformations 
and their dependence on substrate stiffness. More realistic 
planar or thin film geometries will be discussed in a future 
publication. Force balance yields d x a xx + d z a xz = 0. In- 
tegrating over the thickness of the film, using a xz {x, z = 
h) = and a xx (x,z = 0) = Eu x {x, 0), and letting 
a = \ Jo dza xx (x, z), we obtain d x a = Eu x (x,0)/h. In 
the limit h << L, we neglect all z dependence and as- 
sume that the only component of the displacement field 
is u x (x,0) = u(x). Combining then the expression for the 
mean stress, a = Bd x u + (oAp with the boundary condi- 
tion we obtain 



= X 2 



dx 2 



CoAp 



(6) 



where A = ^/Bh/E is a length scale controlled by the 
interplay of cell and substrate stiffness. The solution of 
this equation with boundary conditions a{x = 0) — a{x = 
L) = is 



a(x) = (Ap 1 - 



cosh[(£ - 2x)/2\}\ 
cosh (L/2X) J 



The deformation field is then given by 

(Ap\ sinh[(L- 2ac)/2A] 



u(x) 



B cosh(L/2A) 



(7) 



(8) 



A finite activity Ap generates stresses and deformations 
in the cell, as shown in the top frame of Fig. [l] In an 
isotropic gel, both the stress and the displacement profiles 
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Fig. 2: The stress a(x)/(A[i (dashed line) and displacement 
u{x)B /(Afi (solid line) profiles of a cell on a substrate a con- 
stant stiffness gradient, described by E(x) = Eqx/L are shown 
as functions of the position x inside the cell for X/L = 0.25. 
The profiles are asymmetric and the stress is localized near 
x = L where the stiffness is largest. 



are symmetric about the cell's mid point and the cell is 
uniformly contracted. The deformation is localized near 
the cell's boundaries. The length scale A determined by 
the ratio of cell to substrate stiffness controls the pen- 
etration of the deformation to the interior of the cell. If 
A ~ L, corresponding to a substrate rigidity El ~ Bh/L 2 , 
the active stresses and deformation extend over the entire 
cell. For a cell layer of length 10 /im, thickness 1 fim and 
elastic modulus B ~ 100 kPa, the substrate rigidity pa- 
rameter El can be estimated to be ~ 1 kPa/fim. The 
total deformation A£ — u(0) — u(L) grows with activity 
and is shown in Fig. [I] (bottom frame) as a function of 
X/L ~ \j\jE. The contraction decreases with increasing 
substrate stiffness and saturates to a finite value for soft 
substrates. 

It is also interesting to consider a substrate of varying 
stiffness, as such substrates can be realized in experiments. 
We consider a constant stiffness gradient, corresponding to 
E(x) — E x/L. In this case Eq. (|6| becomes 



XL 



d^a 

dx 2 



1 da 

x dx 



Co A/i 



(9) 



A closed solution can be obtained in terms of hyper- 
geometric functions. The corresponding stress and dis- 
placement profiles are now asymmetric and are shown in 
Fig. [2] The stress is largest in the region of stiffest sub- 
strate, with a correspondingly smaller cell deformation. 
In other words, the largest cell deformation is obtained 
in the boundary region where the substrate is softest. In 
real cells the region where the substrate is softer and the 
resulting stresses in the cell are smaller may correspond 
to region of reduced focal adhesions. Hence the gradi- 
ent stiffness may yield a gradient in the strength of cell- 
substrate adhesion, providing a possible driving force for 
durotaxis, the tendency of cells to move from softer to 



Polarized cell. — We now consider the case of a po- 
larized cell, described by the full free energy /. The cell is 
modeled again as a thin film of length L in the quasi-ld ge- 
ometry described earlier. We are interested in steady state 
configurations. In the chosen geometry these are given by 
the solutions of the equations 
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where P = p(x)x and we have let w' = w and /3/(LT) = 
C3A/1. In the following we scale lengths with the cell's 
length L and stresses with the cell's compressional modu- 

we can eliminate u 
a/B and 



lus B. By combining Eqs. (10a I- (10c 
and rewrite them as coupled equations for a 
p as 

A 2 

a = —^a" + u + v a p 2 + tip' 
L z 



(11a) 



(vp + 2v a w) pp' = Kp" + wa' - [a + bp 2 ) p (lib) 

where the prime denotes a derivative with respect to x/L, 
vo,a,/3 = £o, a ,pA.fi/B, to = 2w/BL, a = a/B, b = b/B, 



and K = K/(BL 2 ) 



to. Thermodynamic stability re- 
there could be 
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stiffcr regions \5p9 30 . 



quires K > 0. As discussed in Ref. 
possible active contributions to the coupling to, which at 
high activity leads to an alternating polarity pattern in 
the gel. Here we restrict ourselves to K > 0. 

In the absence of activity (Ap = 0) Eqs. (11a) and (lib) 
have two homogeneous solutions that satisfy the boundary 
condition a(0) = a{L) = 0, corresponding to an isotropic 
state for a > 0, with p{x) — u{x) — and to a polarized 
state for a < 0, with p(x) = Pa = <J—a/b and u(x) = 0. 
In both cases a(x) = 0. 

For finite activity (A/i 7^ 0), we find two qualitatively 
different solutions, depending on the boundary conditions 
used for the polarization. When Eqs. ( lla[ ) and (lib) are 
solved with boundary condition p(0) = p(L) = 0, consis- 
tent with an isotropic state in the limit A/i = 0, the stress 
is an even function of shown in the top frame of Fig. [3] 
It exhibits a maximum at x — L/2 and is symmetric about 
the mid point of the cell. Both the displacement and the 
polarization vanish at x — L/2 and are odd functions of x 
about this point. For a < we solve the nonlinear equa- 
tions with boundary condition p(0) = p{L) = ^—a/b, 
consistent with a polarized state in the limit A/i = 0. In 
this case the stress, deformation and polarization profiles 
are all asymmetric, as shown in the bottom frame of Fig. [3] 
The sign of the anisotropy is controlled by the sign of the 
polar coupling w. The figure displays the case w > 0, cor- 
responding to filament convection towards the direction of 
positive polarization. 

To quantify the different properties of these two states, 
we define an excess mean polarization averaged over the 
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Fig. 3: Stress a(x)/B (dashed line), deformation field u(x)/L 
(solid line), and polarization Sp(x) = p(x) — po (dotted line) 
profiles obtained by numerical solution of Eqs. ( |lla[ ) and ( |llb[ ) 
for two sets of boundary conditions on the polarization: p(0) = 
p(L) — (top frame) and p(0) = p(L) = 1 (bottom frame). 
Both plots are for X/L — 0.25, w = 4, vq = v a = vp = 1, 
2 = 6=1, # = 1. 



cell as (Sp) = J Q L j^\p(x) — p ]. The excess polariza- 
tion (Sp) is zero for the symmetric polarization profiles 
obtained with the boundary condition p(0) = p(L) = 0, 
whereas (Sp) obtained for the boundary condition p(0) = 
p(L) — \J ~ajb is a non-monotonic function of substrate 
stiffness, as shown in Fig. [4] for three values of activity. 
The excess polarization is largest at a characteristic sub- 
strate stiffness, suggesting that enhancement of stress fiber 
and resulting cell polarization may be obtained for an op- 
timal substrate rigidity, as reported in 



22 . The excess 



polarization (Sp) vanishes in the absence of activity and 
its maximum value increases with activity. 

We have presented a minimal continuum model of the 
interaction of a cell adhering to an elastic substrate. The 
cell is described as an active elastic gel and the coupling to 
the substrate enters as a boundary condition. The model 
shows that the interplay of substrate coupling and activity 
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Fig. 4: Excess mean polarization (Sp) as a function of L/X ~ 
yE obtained from averaging the numerical solutions of Eqs. 



(llal and (fib I for three different values of activity u = vq = 
Va = '■ v = 0.5 (dashed line), v = 1.0 (dotted line) and 
v = 1.5 (solid line). The plots are for w = 4, a — b — 1 and 
K = 1. 



yields contractile stresses and deformation in the cell and 
can enhance polarization, breaking the front/rear symme- 
try of the cell. The model provides a simple, yet powerful 
continuum formulation for the description of cell-substrate 
interactions and can be extended in various directions by 
considering more realistic two-dimensional cell geometries 
and anisotropic and deformable substrates. The possibil- 
ity of cell migration will also be incorporated in future 
work. Finally, the continuum model can be used to de- 
scribe the interaction of confluent layers of epithelial cells 
with substrates. In this case a direct comparison with re- 
cent experiments that have imaged the stress distribution 



in migrating cell layers 31 may be possible 



Acknowledgements This work was supported by 
the National Science Foundation through awards DMR- 
0806511 and NSF-DMR-1004789. We thank Yaouen Fily 
and Silke Henkes for important feedback on cell-substrate 
elasticity and Tannie Liverpool for illuminating discus- 
sions on active systems in general. 



REFERENCES 

[1] Discher D.E., Janmey P. A. and Wang Y., Science, 

310 (2005) 1139. 
[2] Engler A.J., Griffin M.A., Sen S., Bonnemann 

C.G., Sweeney H.L. and Discher D.E., J. Cell Biol, 

166 (2004) 877. 



p-5 



S. Banerjeel M.C. Marchettil,2 



Engler A.J., Sen S., Sweeney H.L. and Discher Zaman M.H., Butler J. P., Weitz D.A., Fredberg 

D.E., Cell, 126 (2006) 677. J.J. and Trepat X., Nature Maternal, 10 (2011) 469. 

Reinhart-King C.A., Methods in enzymology, 443 
(2008) 45. 

Lo CM., Wang H.B., Dembo M. and Wang Y., Bio- 
physical Journal, 79 (2000) 144. 

Yeung T., Georges P.C., Flanagan L.A., Marg B., 
Ortiz M., Funaki M., Zahir N., Ming W., Weaver 
V. and Janmey P. A., Cell motility and the cytoskeleton, 
60 (2005) 24. 

Guo W., Frey M.T., Burnham N.A. and Wang Y., 
Biophys. J., 90 (2006) 2213. 

Barnhart E.L., Lee K.C., Keren K., Mogilner A. 
and Theriot J. A., PLoS Biology, 9 (2011) el001059. 
Kruse K., Joanny J.-F., Julicher F., Prost J. and 
Sekimoto K., Phys. Rev. Lett, 92 (2004) 78101. 
Kruse K., Joanny J.-F., Julicher F., Prost J. and 
Sekimoto K., Eur. Phys. J. E: Soft Matter and Biological 
Physics, 16 (2005) 5. 

Julicher F., Kruse K., Prost J. and Joanny J.-F., 
Phys. Rep., 449 (2007) 3. 

Voituriez R., Joanny J.-F. and Prost J., Europhys. 
Lett., 70 (2005) 118102. 

Giomi L., Marchetti M.C. and Liverpool T.B., Phys. 
Rev. Lett, 101 (2008) 198101. 

Mizuno D., Tardin C., Schmidt C.F. and Mackin- 
tosh F.C., Science, 315 (2007) 370. 
Mackintosh F.C. and Levine A. J., Phys. Rev. Lett, 
100 (2008) 018104. 

Liverpool T.B., Marchetti M.C, Joanny J.-F. and 
Prost J., EurPhys. Lett, 85 (2009) 18007. 
Gunther S. and Kruse K., New Journal of Physics , 9 
(2007) 417. 

Banerjee S. and Marchetti M.C, Soft Matter, 7 
(2011) 463. 

Yoshinaga N., Joanny J.-F., Prost J. and Marcq P., 
Phys. Rev. Lett., 105 (2010) 238103. 
Trepat X., Wasserman M.R., Angelini T.E., Mil- 
let E., Weitz D.A., Butler J. P. and Fredberg J. J., 
Nature Physics, 5 (2009) 426. 

Lima J.I., Sabouri-Ghomia M., Machacek M., Wa- 
terman CM. and Danuser C, Exp. Cell Res., 316 
(2010) 2027. 

Zemel A., Rehfeldt F., Brown A. E.X., Discher 
D.E. and Safran S.A., Nature Physics, 6 (2010) 468. 
Liverpool T.B. and Marchetti M.C, Phys. Rev. Lett., 
90 (2003) 138102. 

Aditi Simha R. and Ramaswamy S., Phys. Rev. Lett, 
89 (2002) 058101. 

Hatwalne Y., Ramaswamy S., Rao M. and Simha 
R.A., Phys. Rev. Lett, 92 (2004) 118101. 
Ahmadi A., Marchetti M.C. and Liverpool T.B., 
Phys. Rev. E, 74 (2006) 061913. 

Kruse K., Joanny J.-F., Julicher F. and Prost J., 
Phys. Biol, 3 (2006) 130. 

Murray J.D. and Oster G.F., Mathematical Medicine 
and Biology, 1 (1984) 51. 

Wong J.Y., Velasco A., Rajagopalan P. and Pham 
Q., Langmuir, 19 (2003) 1908. 

We thank the referee for helping us clarify this point. 
Tambe D.T., Hardin C.C, Angelini T.E., Rajen- 
dran K., Park C.Y., Serra-Picamal X., Zhou E.H., 

p-6 



